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Abstract 

To study the male and female lineages of East Asian and European humans, we have sequenced 25 short tandem repeat markers on 
453 Y-chromosomes and collected sequences of 72 complete mitochondrial genomes to construct independent phylogenetic trees 
for male and female lineages. The results indicate that East Asian individuals fall into two clades, one that includes East Asian 
individuals only and a second that contains East Asian and European individuals. Surprisingly, the European individuals did not 
form an independent clade, but branched within in the East Asians. We then estimated the divergence time of the root of the 
European clade as -41 ,000 years ago. These data indicate that, contrary to traditional views, Europeans diverged from East Asians 
around that time. We also address the origin of the Ainu lineage in northern Japan. 
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Introduction 

The study of human evolution has benefited from the devel- 
opment of genome-wide approaches to obtaining molecular 
data. However, most genes in the genome do not evolve fast 
enough for the study of human evolution at the level of direct 
sequence comparisons, for which reason short tandem repeat 
(STR) markers on the Y chromosome (Y-STR) and mitochon- 
drial DNA sequence (mtDNA), both of which evolve faster 
than autosomal loci, are useful in human evolution studies. 
Those two types of sequences also permit the independent 
study of male and female lineage divergences. Y-STR and 
mtDNA data have thus contributed to elucidating many 
aspects of human evolution (Cann et al. 1987; Excoffier 
et al. 1987; Bowcock et al. 1991; Horai 1991; Excoffier 
et al. 1992; Cavalli-Sforza et al. 1994; Hammer and Horai 



1995; Zerjal et al. 1997; Bergen et al. 1999; Shinka et al. 
1999; Tanaka et al. 2004; Katoh et al. 2005a; Derenko 
et al. 2007). In addition, while Ingman et al. (2000) and 
Maca-Meyer et al. (2001) have shown that mtDNA sequences 
are still appropriate for the studies of global human diversity, 
Katoh et al. (2005b) demonstrated the usefulness of Y-STR for 
the evolutionary study of East Asian ethnic groups. 

However, those studies have lacked resolution for individ- 
ual lineages and have not provided coherent pictures of male 
and female lineage splittings. Indeed, a major issue in human 
evolution is that the study of the male lineage delivers results 
that differ from those obtained through the study of the 
female lineage (Underhill and Kivisild 2007). Therefore, we 
focused in the present study on examining human evolution 
at the level of males and females separately. Our interest 
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focused on the divergence of three major populations — 
Africans, Europeans, and East Asians — and East Asians and 
Europeans in particular, because they lived together for a 
while after coming out of Africa and likely interbred, perhaps 
similar to the hybridizations reported between Neanderthals 
and modern humans (Noonan 2006). After the divergence 
between Mongoloid (East Asian) and Caucasoid (European) 
people -55,000 years ago (Nei and Roychoudhury 1 993), they 
might have still interbred while en route to settlement in their 
present localities. 

To investigate this, we sampled a number of European and 
East Asian individuals to obtain information on the male and 
female lineages, by constructing phylogenetic trees based on 
Y-STR and mtDNA sequences, respectively. These two trees 
should provide new information on the divergence of the East 
Asian and European lineages. We included in our study data 
from the Ainu people in the northmost island of the Japanese 
archipelago, because their ancestor and evolutionary history 
are still obscure. Recently, Jinam et al. (2012) suggested that 
the Ainu are closer to the people in Okinawa from the south- 
most archipelago in Japan than they are to any other East 
Asian people. Our investigation also addresses Ainu origin. 

Materials and Methods 

Population Samples 

To study male lineages, we sequenced and analyzed Y-STR for 
453 unrelated males, which were classified into 242 Japanese, 
85 Korean, 77 Khalkh Mongolian, and 49 European bins. The 
Korean samples were collected in Seoul, while the Khalkh 
Mongolian samples were obtained in Ulaanbaatar, Mongolia 
(Katoh et al. 2002, 2005a, 2005b). The Khalkh are the major 
population in Mongolia, representing ~ 80% of the total pop- 
ulation, the language spoken in Mongolia has its origins in 
the Khalkh language. The European samples were collected 
from the Coriell Cell Repositories (Human variation panel, 
Caucasoid, HD100CAU). Informed consents were obtained 
from donors at each laboratory or institution. For the genetic 
study of these samples, we obtained permission from the 
ethics committee of Tokai University, Japan. 

Y Chromosome Binary Markers 

We typed 453 male individuals for a total of 1 4 binary markers 
on the Y chromosome. These markers were M130, YAP, 
M15, M57, SRY-4064, M89, M9, M128, Tat, M175, M119, 
SRY+465, M122, and M45. M130 was typed by polymerase 
chain reaction (PCR) amplification using the primers described 
in Bergen et al. (1999), followed by Bsll digestion. YAP, SRY- 
4064, Tat, and SRY+465 were typed according to Hammer 
and Horai (1995), Bravi et al. (2000), Zerjal et al. (1997), and 
Shinka et al. (1999), respectively. 

For the indel markers, M15, M57, M128, and M175, we 
performed PCR amplification with the primers described in 



Underhill et al. (2001), labeling the 5' end of the forward 
primer with a fluorescent marker. For M89, M9, M122, 
M119, and M45, we carried out allele-specific PCR assays 
according to the primer designs of Su et al. (1999). These 
markers consisted of two forward primers, which defined 
SNP alleles, and one reverse primer. We labeled the 5 r end 
of each forward primer with a different fluorescent dye (FAM 
and HEX). PCR reactions were performed in a 10-jil reaction 
volume, containing 1x AmpliTaq PCR buffer with 1.5mM 
MgCI 2 (Applied Biosystems), 200 jiM each of dNTP, 0.25 U 
AmpliTaq polymerase (Applied Biosystems), 1 pmol of each 
primer, and 6ng of genomic DNA. The cycle conditions 
were as follows: 94°C for 2min; followed by 30 cycles of 
94 °C for 30 s, 54-62 °C for 30 s, and 72 °C for 50 s; followed 
by a 5-min extension at 72 °C. The amplified products were 
detected using the ABI Prism 3700 DNA analyzer. 

Y-STR Genotyping 

We used 25 STR markers on the Y chromosome: DYS388, 
DYS389I, DYS389b, DYS390, DYS391, DYS392, DYS393, 
Y-GATA-A7.1 (DYS460), Y-GATA-A7.2 (DYS461), Y-GATA- 
A10, Y-GATA-C4, Y-GATA-H4, DYS597, DYS598, DYS599, 
DYS600, DYS601, DYS602, DYS603, DYS604, DYS605, 
DYS606, DYS607, DYS608, and DYS609. We labeled the 
5' end of each forward primer with a fluorescent dye (FAM 
or HEX), and a PCR reaction for each marker was carried 
out independently. PCR reactions for DYS388, DYS389I, 
DYS389b, DYS390, DYS391, DYS392, and DYS393 were per- 
formed according to published procedures with minor modi- 
fications (Kayser et al. 1 997). PCR reactions for the remaining 
18 markers were carried out in a total volume of 10 con- 
taining 1 0 mM Tris-HCI, pH 8.0, 50 mM KCI, 2.0 mM MgCI 2 , 
0.01% gelatin, 200 jiM each of dNTP, 0.3 U AmpliTaq poly- 
merase (Applied Biosystems), 2 pmol of each primer, and 6 ng 
of genomic DNA. These markers were amplified by the 
following cycling conditions: pre-PCR of 5min at 96°C, 
1 min at 57 °C, and 1 min at 72 °C; followed by 30 cycles of 
96 °C for 45 s, 57 °C for 45 s, and 72 °C for 1 min; with a final 
extension at 72 °C for 5 min. 

Amplicons were separated using the ABI Prism 3700 DNA 
analyzer, and allele sizes were calibrated with the GeneScan 
500 ROX size standard. Fragment sizes were determined using 
GeneScan Analysis 3.5 and Genotyper 3.5NT software 
(Applied Biosystems). The number of repeat units in each frag- 
ment was determined by sequencing reference DNA samples 
with different alleles. The alleles were designated with repeat 
numbers. The allele length for DYS389b was obtained by sub- 
tracting the allele length of DYS389I from that of DYS389II 
(Qamar et al. 2002; Zerjal et al. 2002). 

Data Analysis 

Y-chromosome binary haplogroups for the four human pop- 
ulations mentioned earlier were defined by the analysis of all 
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Fig. 1 . — Phylogenetic tree of Yap haplotypes for Europeans and East Asians. The tree was constructed by using the Y-STR data, and the Yap haplotypes 
were placed at the tips of the tree. In the tree, the haplotypes are classified into Yap-A and Yap-B first, and Yap-B is further divided into Yap-B 1 and Yap-B2. 
Yap-B 1 includes East Asian individuals only, while Yap-B2 contains European and East Asian individuals together. 



14 binary markers. The nomenclature of haplogroups fol- 
lowed that of the Y chromosome consortium (Slatkin 1995) 
and Karafet et al. (2008). Then, a Y-STR haplotype tree for the 
453 male individuals and a mtDNA tree for the 72 individuals 
were constructed by using the neighbor-joining (NJ) method 
(Saitou and Nei 1 987) in the MEGA 5.0 (Tamura et al. 201 1). 
The distance computed for the Y-STR tree was the R S j dis- 
tance (Slatkin 1995) between a pair of Y-STR haplotypes. For 
haplotypes x and y, we computed the distance (d) between 
them by using the formula 

d = J2(x l -y l ) 2 /n, 

/=1 

where n is the number of loci, x, and y- are the number of 
repeats of the /th locus of the haplotypes x and y, respectively 
(fig. D. 

To study the female lineages we collected the complete 
mtDNA genomes (16,750 bp each) from the International 
Nucleotide Sequence Data Collaboration (DDBJ/ EMBL/ 
GenBank) to study the female lineage of the European and 
East Asian. The mtDNA genome data include 32 Japanese, 
1 1 Korean, 5 Mongolian, 23 European, and 1 African samples 
(table 1). For the mtDNA tree the evolutionary distance 
between a pair of individuals was obtained by the Kimura 
two-parameter method (Kimura 1980). The bootstrap test 
(Felsenstein 1985) was performed for 1,000 replicates for 
each tree (fig. 2). The African sample was used as the 
reference in the phylogenetic trees. 



Results 

For the male lineages, we used the Y-STR markers for 453 
individuals covering Japanese, Korean, Mongolian (Khalkh), 
American, and European people. We determined the evolu- 
tionary distances (R S j) among them, and constructed a phylo- 
genetic tree using the NJ method (Saitou and Nei 1987), as 
shown in figure 1 . The Y-STR tree revealed that the male an- 
cestral lineage contained two clades (Yap-A and Yap-B). While 
Yap-A clade includes the East Asian individuals only, Yap-B 
clade contains the East Asian and European individuals 
together. Surprisingly, the European males never formed an 
independent clade. Instead, they formed separate clades 
within Yap-B. 

We then constructed a phylogenetic tree (Mt tree) for 
the 72 complete mtDNAs including the four ethnic groups, 
as shown in figure 2. We used Kimura's two parameter 
method (Kimura 1980) for computing evolutionary distances 
among them and the NJ method for the tree construction. The 
Mt tree revealed two female descendant clades (Mt A and Mt 
B). Mt B consisted of the East Asian females only, while Mt A 
contained the East Asian and European females together. 
As in the case of males, the European females did not form 
an independent clade, but comprised several groups within 
the Mt A tree. As the Mt A cluster includes roughly as many 
European individuals as East Asian individuals, it is not clear 
which of them is ancestral to the other. The node marked with 
the blue circle in figure 2 suggests that the East Asians are 
ancestral to the Europeans. The bootstrap value of the node 
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Table 1 

List of Accession Number for Mitochondrial Genomes Analyzed in This Study 



Group Name 



Accession Number 



Japanese 32 AP008912, AP010835, AP010837, AP008691, AP010664, AP008279, AP008552, AP008751, AP010834, AP008837, AP010832, 

AP008791, AP008276, AP008270, AP008280, AP008274, AP008705, AP008301, AP008283, AP008777, AP010831, AP008267, 
AP008272, AF346990, AP010663, AP010676, AP008768, AP008302, AP008303, AP010830, AP008706, AP008627 

Caucasian 23 EU600323; France Druze, JQ343921; France, GQ200588; USA, GU294854; Canada, HQ399469; USA, JQ048704; Portugal, 
KC618506; England: Devon, NC 012920; Great Britain, FJ216960; Europe, HQ729918; France: Brittany, HQ610202; USA, 
GU323604; Ireland, KC469897; England Worcestershire, KC862290; England: Liverpool, GQ129173; Spain: Southern Spain, 
AY275537; Spain: Canary Islands, AY275530; Spain: Galicia, HQ025914; USA, HQ651683; Portugal, HQ651702; Portugal, 
HQ651697; Portugal, HQ585390; USA, HQ651686; Portugal 

Korean 11 AF346993, EF1 53824, EF1 53823, EF397561, FJ951594, EF1 53822, EF1 53821, FJ951589, FJ951593, FJ951590, FJ951592 

Mongolian 5 EU007891, EU007892, AY255146, EU007893, EU007890 

African 1 AF346995 



is 99%. Therefore, both male and female lineages suggest 
that Europeans diverged from within East Asian ancestors or 
that they interbred with East Asian individuals up to a certain 
divergence time. 

Our next question was thus to estimate the divergence time 
of the European clade within the East Asian lineage, for males 
and females separately. To address that question, we com- 
puted the evolutionary distance (R SJ ) between every pair of 
the male individuals to construct the Y-STR tree in figure 1 . 
The R S j value between the bottom and root of A and B clades 
in the tree was 16.91, while that between the bottom and 
root of the European male clade in the B clade was 12.31. 
Based on a divergence of East Asians from Africans of 55,000 
years ago (Nei and Roychoudhury 1974, 1993), and assuming 
that R SJ is proportional to time, we can estimate the evolu- 
tionary rate of Y-STR by using the Y-STR tree. In the tree the 
R S j value between the deepest root and the bottom is 16.91, 
and that between the common ancestor of Europeans males 
and the bottom is 12.31 . The rate is thus estimated as 16.91/ 
55,000 = 3.07 x 10~ 4 per repeat per year, which leads us to 
the conclusion that the divergence time of the European 
males is -40,100 years ago. 

In the case of the Mt tree, we first transformed it into tree 
in which the lengths of a pair of branch lengths from the com- 
mon node were equal to one another, as in the UPGMA tree 
(Sokal 1958), because we now dealt with the evolutionary 
time rather than distance (Kumada 1993). In the transformed 
tree, the topology was kept unchanged. We call the trans- 
formed tree the evolutionary time (ET) tree (fig. 3), because 
the tree reflects the evolutionary time rather than the distance, 
we call the tree showing evolutionary distance the evolution- 
ary distance (ED) tree. 

In the present ET tree, the distance between the bottom 
and root in any route was 0.00132 per nucleotide site, while 
that between the bottom and the root (marked with the 
closed circle) of the European male clade was 0.000996. 
Therefore, assuming again that the two major clades in the 



Mt tree diverged from the African -55,000 years ago, we 
estimated the evolutionary rate of mtDNA as 2.4 x 10~ 8 per 
site per year. Our estimate was comparable with that of 
Ingman et al. (2000). It is noted that our estimate for the 
human mitochondrial genome is more than ten times faster 
than that of the human nuclear neutral sequences (Fukami 
Kobayashi 2005). Then, by applying the evolutionary rate of 
mtDNA to the distance between the bottom and root in the 
European clade, we estimated the divergence time of the 
European females, as -41,500 years ago. 

Thus, both the Y-STR and Mt trees are mutually consistent 
with respect to the tree shape and divergence time of 
European individuals relative to East Asians, despite being 
based on independent data and distance measures. Our re- 
sults suggest that the European people settled down in their 
territories -41,000 year ago, and have developed their own 
cultures and languages since then. 

On the other hand, the East Asians were classified into two 
clusters; one is Type 1 East Asians denoted as Yap-A and Mt B 
clusters in the male and female trees, respectively, and the 
other is Type 2 East Asians denoted as Yap-B and Mt A clusters 
in the male and female trees, respectively. While Type 1 in- 
cludes East Asian individuals only, Type 2 contains East Asians 
and European individuals together. The general view of 
the East Asian and European divergences is summarized in 
figure 4. 

Discussion 

Both the Y-STR and mtDNA trees consistently show that 
Europeans diverged from East Asian ancestors -41,000 
years ago. Population genetic theory indicates that 41,000 
years, or about 2,000 generations, are long enough to ac- 
cumulate SNPs in the same loci in each lineage (Kimura 
1 983; Nei 1 987) to account for the present genetic and phe- 
notypic differences between the East Asians and Europeans, 
but too short to acquire independent loci between them. 
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gi|316992556|gb|HQ651686.1 1 Portugal 
gi|312231545|gb|HQ585390.1| ethnicity:Caucasian; originjocality: USA 
gi|31 699271 0|gb|HQ651 697.1 1 Portugal 
gi|316992780|gb|HQ651702.1| Portugal 
gi|316992514|gb|HQ651683.1| Portugal 

gi|303325296|gb|HQ025914.1| ethnicity:Caucasian; originjocality: USA 
|33465939|gb|AY275530.1| Spain: Galicia 

gi|395764494|gb|AY275537.2| Spain: Canary Islands 
gi|237652923|gb|GQ1 29173.1 1 Spain: Southern Spain 

gi|478343666|gb|KC862290.1 1 ethnicity:English; originjocality: England: Liverpool 
gi|154101200|gb|EU007891.1 1 Mongolian 

gi|442565727|gb|KC469897.1| ethnicity:English; origin_locality:England Worcestershire 
gi|283049431 |gb|GU323604.1 1 ethnicity:Caucasian; origin_locality:lreland 
gi|31 2231 601 |gb|HQ61 0202.1 1 ethnicity:Caucasian; origin_locality:USA 
■gi|316926349|gb|HQ729918.1| French; origin_locality:France: (Coatreven) Brittany 
gi|206575353|gb|FJ21 6960.1 1 ethnicity:Caucasian; originjocality: Europe 
gi|251831106|ref|NC 012920.1 1 United Kingdom: Great Britain 
i|451686407|gb|KC618506.1| ethnicity:English; originJocality:England: Devon 
gi|13272976|gb|AF346993.1| Korean 
gi|61288301 |dbj|AP008912.1 1 Japan Tokyo 
gi|220897773|dbj|AP01 0835.1 1 Japan: Chiba 

gi|220897801|dbj|AP010837.1 1 Japan: Tokyo 
gi|61287925|dbj|AP008691.1 1 Japanese male Aichi 
gi|61287256|dbj|AP008279.l| Japan: Tokyo 
gi 1 1 71 452437|d bj | AP01 0664.1 1 Japa n f ema le 
gi|154101214|gb|EU007892.1 1 Mongolian 

gi|359295022|gb|JQ048704.1| Portuguese; originjocality: Portugal 

gi|309400117|gb|HQ399469.1| ethnicity:Caucasian; originjocality: USA 

gi|61288038|dbj|AP008751.1| Japan Aichi 

gi|61287696|dbj|AP008552.1 1 Japan Aichi 
gi|61288172|dbj|AP008837.l| Japan: Tokyo 
gi|220897759|dbj|AP01 0834.1 1 Japan: Kanagawa 

gi|281 1 901 22|gb|GU294854.1 1 ethnicity:Caucasian; originjocality: Canada 
gi|238627114|gb|GQ200588.1| ethnicityxaucasian; originjocality: USA 

gi|371491686|gb|JQ343921.1| ethnicity:French; originjocality: France 
gi|187960893|gb|EU600323.1| France Druze 

gi|220897731 |dbj|AP01 0832.1 1 Japan: Ehime 



MtA 



Egi|61288099|dbj|AP008791.l| Japan Aichi 
gi|119035558|gb|EF153824.1| South Korea: Seoul 
— gi |61 287251 |dbj| AP008276.1 1 Japan Tokyo 
gi 1 1 24388042 |gb| E F397561 .1 1 Korean 
gi|1 1 9035544|gb|EF1 53823.1 1 Korean 
gi |1 541 01 1 86|gb| EU007890.1 1 Mongolian 



gi|32891 271 |gb|AY2551 46.1 1 Mongolian from Inner Mongolia 
gi|61 287241 |dbj|AP008270.1 1 Japan Tokyo 



gi|61287249|dbj|AP008274.l| Japan Tokyo 



gi|61287257|dbj|AP008280.l| Japan Tokyo 

gi|229503574|gb|FJ951594.1| South Korea: Seoul 



gi|61287969|dbj|AP008705.1 1 Japan Aichi 



88 | 



gi|119035516|gb|EF153821.l| South Korea: Seoul 
— gi|1 1 9035530|gb|EF1 53822.1 1 South Korea: Seoul 

gi|61287294|dbj|AP008301.l| Japan Tokyo 

gi|61288079|dbj|AP008777.1| Japan Aichi 



gi|61287263|dbj|AP008283.1| Japan Tokyo 
- gi|229503504|gb|FJ951589.1 1 Korean 

gi|154101228|gb|EU007893.1| Mongolian 

gi|61287237|dbj|AP008267.1| Japan Tokyo 

gi|220897717|dbj|AP010831.1| Japan: Okayama 

gi|61287245|dbj|AP008272.1 1 Japan Tokyo 



MtB 



90 J 



- gi|61 28781 9|dbj|AP008627.1 1 Japan Aichi 

- gi|61287970|dbj|AP008706.1| Japan Aichi 
gi|229503546|gb|FJ951592.1| Korean 




gi|171452423|dbj|AP01 0663.1 1 Japan female 
gi|13272934|gb|AF346990.1 1 Japanese individual 

gi|171452605|dbj|AP01 0676.1| Japan female 



gi|229503560|gb|FJ951 593.1 1 Korean 



Red: European 
Blue: East Asian 



gi|220897703|dbj|AP01 0830.1 1 Japan: Aichi 

gi|61287298|dbj|AP008303.1| Japan Tokyo 
gi|61287297|dbj|AP008302.l| Japan Tokyo 
gi|229503518|gb|FJ951590.1| Korean 
gi|61288064|dbj|AP008768.1| Japan Aichi 

gi|13273004|gb|AF346995.1 1 African (Mandenka) 



Fig. 2. — Phylogenetic tree of the complete mitochondrial genomes for the four ethnic groups. It is divided into MtA and MtB. MtA contains East Asian 
and European individuals together while MtB includes East Asian individuals only. 
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31 6992556|gb|HQ651 686. 1 1 Homo sapiens isolate P6 mitochondrion complete genome 
31 2231 545|gb|HQ585390. 1 1 Homo sapiens haplogroup U6a3 mitochondrion complete genome 
|31 699271 0|gb|HQ651 697.1 1 Homo sapiens isolate P17 mitochondrion complete genome 
31 6992780|gb|HQ651 702. 1 1 Homo sapiens isolate P22 mitochondrion complete genome 
|31 699251 4|gb|HQ651 683.1 1 Homo sapiens isolate P3 mitochondrion complete genome 
|303325296|gb|HQ025914.1 1 Homo sapiens haplogroup U6a7 mitochondrion complete genome 
|33465939|gb|AY275530.1| Homo sapiens isolate 1766 mitochondrion complete genome 
|395764494|gb|AY275537.2| Homo sapiens isolate 183 mitochondrion complete genome 
237652923|gb|GQ1 291 73. 1 1 Homo sapiens isolate 42 U5b3d (Tor289) mitochondrion complete genome 
|478343666|gb|KC862290.1| Homo sapiens haplogroup U5a1a2a1 mitochondrion complete genome 
154101200|gb|EU007891 .1 1 Homo sapiens isolate 220 mitochondrion complete genome 
|442565727|gb|KC469897.1| Homo sapiens haplogroup K2b1 mitochondrion complete genome 
283049431 |gb|GU323604. 1 1 Homo sapiens haplogroup K1c2 mitochondrion complete genome 
|312231601 |gb|HQ61 0202.1 1 Homo sapiens haplogroup K1a3a mitochondrion complete genome 
31 6926349|gb|HQ72991 8. 1 1 Homo sapiens haplogroup V mitochondrion complete genome 
|206575353|gb|FJ21 6960.1 1 Homo sapiens haplogroup H1 mitochondrion complete genome 
251 831 1 06|ref|NC 012920.1 1 Homo sapiens mitochondrion complete genome 
451 686407|gb|KC61 8506.1 1 Homo sapiens haplogroup H33 mitochondrion complete genome 
|13272976|gb|AF346993.1| Homo sapiens mitochondrion complete genome 
61 288301 |dbj| AP00891 2. 1 1 Homo sapiens mitochondrial DNA complete genome isolate: JDsq0096 
|220897773|dbj|AP01 0835.1 1 Homo sapiens mitochondrial DNA complete genome isolate: CB1068 
|220897801 |dbj|AP01 0837.1 1 Homo sapiens mitochondrial DNA complete genome isolate: TK1079 
61 287925|dbj|AP008691 . 1 1 Homo sapiens mitochondrial DNA complete genome isolate: HNsq0215 
61 287256|dbj|AP008279. 1 1 Homo sapiens mitochondrial DNA complete genome isolate: TCsq0020 
|171452437|dbj|AP010664.1| Homo sapiens mitochondrial DNA complete genome isolate: SCsq0004 
|154101214|gb|EU007892.1| Homo sapiens isolate 221 mitochondrion complete genome 
|359295022|gb|JQ048704.1| Homo sapiens haplogroup J1c3b mitochondrion complete genome 
3094001 1 7|gb|HQ399469. 1 1 Homo sapiens haplogroup T2b mitochondrion complete genome 
61 288038|dbj|AP008751 . 1 1 Homo sapiens mitochondrial DNA complete genome isolate: NDsq0056 
|61287696|dbj|AP008552.1| Homo sapiens mitochondrial DNA complete genome isolate: ONsq0025 
|61288172|dbj|AP008837.1| Homo sapiens mitochondrial DNA complete genome isolate: JDsq0015 
220897759 |d bj | A P 0 1 0834.1 1 Homo sapiens mitochondrial DNA complete genome isolate: KN221 
281 1 901 22|gb|GU294854. 1 1 Homo sapiens haplogroup 12 mitochondrion complete genome 
2386271 14|gb|GQ200588.1 1 Homo sapiens haplogroup X2 mitochondrion complete genome 
|371491686|gb|JQ343921.1| Homo sapiens haplogroup X2b4 mitochondrion complete genome 
|187960893|gb|EU600323.1| Homo sapiens isolate Druze292 mitochondrion complete genome 
|220897731 |dbj|AP01 0832.1 1 Homo sapiens mitochondrial DNA complete genome isolate: EH1061 
61 288099|dbj|AP008791 . 1 1 Homo sapiens mitochondrial DNA complete genome isolate: NDsq0157 
119035558|gb|EF153824.1| Homo sapiens isolate Korean 97 mitochondrion complete genome 
61287251 |dbj|AP008276.1 1 Homo sapiens mitochondrial DNA complete genome isolate: TCsq0017 
124388042|gb|EF397561.1| Homo sapiens isolate Korean 19 mitochondrion complete genome 
|119035544|gb|EF1 53823.1 1 Homo sapiens isolate Korean 53 mitochondrion complete genome 
154101 186|gb|EU007890.1 1 Homo sapiens isolate 219 mitochondrion complete genome 
32891 271 |gb|AY2551 46. 1 1 Homo sapiens isolate Mg246 mitochondrion complete genome 
61287241 |dbj|AP008270.1 1 Homo sapiens mitochondrial DNA complete genome isolate: TCsqOOH 
|61287249|dbj|AP008274.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0015 
|61287257|dbj|AP008280.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0021 
229503574|gb|FJ951 594. 1 1 Homo sapiens isolate Korean Kor75 mitochondrion complete genome 
|61287969|dbj|AP008705.1| Homo sapiens mitochondrial DNA complete genome isolate: HNsq0231 
|119035516|gb|EF153821.1| Homo sapiens isolate Korean 33 mitochondrion complete genome 
1 19035530|gb|EF1 53822.1 1 Homo sapiens isolate Korean 38 mitochondrion complete genome 
61 287294|dbj|AP008301 . 1 1 Homo sapiens mitochondrial DNA complete genome isolate: TCsq0042 
61 288079|dbj|AP008777. 1 1 Homo sapiens mitochondrial DNA complete genome isolate: NDsq0117 
61 287263|dbj|AP008283. 1 1 Homo sapiens mitochondrial DNA complete genome isolate: TCsq0024 
229503504|gb|FJ951 589. 1 1 Homo sapiens isolate Korean Kor1 mitochondrion complete genome 
1 541 01 228|gb|EU007893. 1 1 Homo sapiens isolate 222 mitochondrion complete genome 
|61287237|dbj|AP008267.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0008 
22089771 7|dbj|AP01 0831 . 1 1 Homo sapiens mitochondrial DNA complete genome isolate: OY1056 
|61287245|dbj|AP008272.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0013 
|171452423|dbj|AP010663.1| Homo sapiens mitochondrial DNA complete genome isolate: SCsq0003 
13272934|gb|AF346990.1| Homo sapiens mitochondrion complete genome 

|171452605|dbj|AP010676.1| Homo sapiens mitochondrial DNA complete genome isolate: SCsq0016 
229503560|gb|FJ951 593. 1 1 Homo sapiens isolate Korean Kor24 mitochondrion complete genome 
|61287819|dbj|AP008627.1| Homo sapiens mitochondrial DNA complete genome isolate: ONsq0142 
|61287970|dbj|AP008706.1| Homo sapiens mitochondrial DNA complete genome isolate: HNsq0232 
229503546|gb|FJ951 592. 1 1 Homo sapiens isolate Korean Kor23 mitochondrion complete genome 
220897703 1 d bj | A P 0 1 0830.1 1 Homo sapiens mitochondrial DNA complete genome isolate: AC1044 
|61287298|dbj|AP008303.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0044 
|61287297|dbj|AP008302.1| Homo sapiens mitochondrial DNA complete genome isolate: TCsq0043 
22950351 8|gb|FJ951 590. 1 1 Homo sapiens isolate Korean Kor4 mitochondrion complete genome 
|61288064|dbj|AP008768.1| Homo sapiens mitochondrial DNA complete genome isolate: NDsq0103 
|13273004|gb|AF346995.1| Homo sapiens mitochondrion complete genome 



Fig. 3. — A phylogenetic tree of the complete mitochondrion genome (16,750 bp) including D-loop regions of 72 humans; 32 Japanese, 1 1 Korean, 
5 Mongolian, 23 European, and 1 African samples. This phylogenetic tree was constructed using the UPGMA method by the MEGA 5.0 version program. 
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Fig. 4. — Phylogenetic tree of the major human populations Africans, 
Type 1 East Asians, Europeans, and Type 2 East Asians. 

Recently, Liu (2012) reported on five genes responsible for 
the facial morphology of European people. The East Asian 
people must have the counterparts that differ at the SNP 
level from those in the European people. As our phylogenetic 
trees demonstrate, the European alleles at the five loci have 
diverged from the ancestral East Asian alleles. 

Our result contrasts with the traditional view that 
Europeans and East Asians simultaneously diverged from 
African ancestors 55,000 years ago. It is noteworthy, how- 
ever, that Shinoda (2007) investigated into the haplogroups of 
mtDNA, and revealed a number of evolutionary haplotype 
lineages. The lineages include L3 (African), N (East Asian), W 
(European), and L3 to N to R (East Asian) and then HV 
(European) among others. Though they did not explain their 
results, their haplotype lineages can now be understood by 
our finding that the Europeans diverged from the East Asians. 
Therefore, the discrepancy between the traditional view and 
ours lies mainly in that the traditional view was based on 
autosomal genes that evolved much slower than Y-STR or 
mtDNA, and could not distinguish the evolutionary lineages 
at the individual level. Note that, as estimated earlier, the evo- 
lutionary rate of mtDNA is 2.4 x 1 0~ 8 per site per year, while 
that of nuclear neutral sequences is 2.0 x 1CT 9 per site per 
year (Fukami Kobayashi 2005). The discrepancy also is due 
to the fact that while we dealt with many male and female 
individuals in our study, the other studies did not. 

By typing the 14 binary markers of the Y-STR sequences 
according to the classification agreed at the Y Chromosome 
Consortium (2002) and elsewhere (Karafet et al. 2008), we 
found that all individuals in the Yap-B1 in figure 1 belonged to 
either haplogroup C or D, while the majority in the Yap- B2 
belonged to haplogroup O. Since Yap-B1 includes mainly 
Japanese and Korean males, in which 91% individuals share 
haplotype 02b, the Korean and Japanese males are definitely 
closest to one another within East Asian humans. We also 
found a high frequency of the 02b haplotype in Manchu 
(Northern China) and Korean-Chinese samples (Katoh et al. 
2005b; Kim et al. 2011). 

The origin of the Ainu people is still an unresolved issue 
(Tajima et al. 2004). On the basis of our results, we propose a 
possible scenario for the origin of the Ainu people, who now 
live in the north-most island of the Japanese archipelago, 
Hokkaido. The Ainu people have European phenotypic 



characters, but they are genetically closer to East Asians 
than to Europeans (Watanabe 1 975). These contradictory fea- 
tures of the Ainu people are puzzling. As shown in figure 4, 
Europeans may have diverged from East Asians -41 ,000 years 
ago, it is possible that hybrid individuals were born before 
the divergence, and some of them looked more like the 
Europeans while possessing a generally East Asian genotype. 
We suggest that the ancestor of the Ainu people was such a 
group of the hybrid individuals. We note that the present Ainu 
people share the mtDNA haplotype not with the Japanese but 
with the European living in Siberia, Russia (Adachi et al. 2009). 
Thus, we furthermore suggest that the ancestor of the Ainu 
originated in northern Eurasia and took a route through 
Siberia and north China before settling in northern regions 
of Japan and nearby places. There is a report that other 
people than the Ainu also lived in the northern regions but 
disappeared (Adachi et al. 2009). As the people in Okinawa 
islands are closest to the Ainu people in the East Asians (Jinam 
et al. 201 2), they might also be descendants from of mixing of 
East Asian and European lineages. 

Conclusion 

Humans are traditionally viewed as falling into three major 
populations, Africans, East Asians, and Europeans, with the 
latter two diverging from the African ancestors -55,000 years 
ago. We sequenced Y-STR makers in Y chromosome and col- 
lected complete mtDNAs for many East Asian and European 
individuals to reexamine that view. Phylogenetic trees of 
Y-STR makers and Mt genes suggest that the Europeans in- 
terbred with East Asians until -41,000 years ago. On the 
other hand, East Asians diverged from their African acestors 
-55,000 years ago. Therefore, we suggest that the European 
and East Asian lineages diverged -41,000 years ago. 
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